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ABSTRACT 

The best-fitting 2-dimensional plane within the 3-dimensional space of spi- 
ral galaxy disk observables (rotational velocity v TO t, central disk surface brightness 
= —2.5 log Io, and disk scale-length h) has been constructed. Applying a three- 
dimensional bisector method of regression analysis to a sample of ^100 spiral galaxy 
disks that span more than four mag arcsec -2 in central disk surface brightness yields 

V mt OC ig-6ff±0.°5fc0.77±0.07 (£>_ band ) and Urot ^ j0-43i0.04jjQ.69i0.07 (i? _ band ). C on- 

trary to popular belief, these results suggest that in the -B-band, the dynamical mass- 
to-light ratio (within 4 disk scale-lengths) is largely independent of surface brightness, 



varying as Iq 



o.ooio.io 



jjO.54io.i4_ Consistent results were obtained when the expanse of 



the analysis was truncated by excluding the low surface brightness galaxies. Previ- 

— 1/2 

ous claims that M/Lb varies with I are shown to be misleading and/or due to 
galaxy selection effects. Not all low-surface-brightness disk galaxies are dark matter 
dominated. The situation is however different in the near-infrared where Lk> oc u 4 

— 1/2 

and M/Lk' is shown to vary as I . Theoretical studies of spiral galaxy disks should 
not assume a constant M/L ratio within any given passband. 

The i?-band dynamical mass-to-light ratio (within 4 disk scale-lengths) has no 
obvious correlation with (B — R) disk colour, while in the if' -band it varies as - 
1.25±0.28(B — R). Combining the present observational data with recent galaxy model 
predictions implies that the logarithm of the stellar-to-dynamical mass ratio is not a 
constant value, but increases as disks become redder, varying as 1.70±0.28(£? — R). 

Key words: cosmology: dark matter - galaxies: fundamental parameters - galaxies: 
kinematics and dynamics - galaxies: photometry - galaxies: spiral 



1 INTRODUCTION 

Given the large observed range in the disk parameters of spi- 
ral galaxies, the tightness (and slope) of the Tully & Fisher 
(1977, hereafter TF) relation between absolute magnitude 
and the logarithm of the rotational velocity is a fundamental 
clue as to how these systems were constructed (Dalcanton, 
Spergel & Summers 1997; Mo, Mao & White 1998; Mo & 
Mao 2000; van den Bosch 2000; Koda, Sofue & Wada 2000). 
Not only does this relationship hold for the High Surface 
Brightness (HSB) spiral galaxies that typify the Hubble se- 
quence, but also for the disk-dominated Low Surface Bright- 
ness (LSB) spiral galaxies (Zwaan et al. 1995). The univer- 
sality of the TF relation for both (HSB) and LSB galaxies 
immediately tells us that for a given rotational velocity one 
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has a constant galaxy luminosity, albeit within the scatter 
of the TF relation. This had been used to infer that, at 
any given constant luminosity, the mass-to-light (M/L) ra- 
tio must vary linearly with increasing disk scale-length h, 
or equivalently as 7 , where Io is the central disk sur- 
face brightness in intensity units (Zwaan et al. 1995). De- 
spite popular misconception, how the M/L ratio may vary 
between galaxies having different luminosities was not ad- 
dressed and is the subject of the present analysis. This paper 
derives the mean dependency of the M/L ratio on the scale- 
length and (central) surface brightness of a galaxy's disk. 

The issue above is tackled by replacing the magnitude 
term in the TF diagram with the central disk intensity and 
scale- length terms (recall, Ldisk = 2irIoh 2 ). This gives three 
fundamental variables: the rotational velocity of the disk 
(frot), the central disk surface brightness (/io = —2.5 log Io 
mag arcsec -2 ), and the disk scale-length (h). Just as the 
Faber- Jackson (1976) relation between absolute magnitude 



© 2001 RAS 



2 Alister W. Graham 



and velocity dispersion for Elliptical galaxies was supple- 
mented through the construction of the 'Fundamental Plane' 
(Djorgovski & Davis 1987; Dressier et al. 1987; Faber et al. 
1987), this study will explore the question of whether or not 
the use of three variables, rather than two are required for 
the treatment of spiral galaxy disks. 

In recent years, numerous authors have addressed this 
question using different galaxy samples and different meth- 
ods of data analysis (e.g., Willick et al. 1997). It should be 
noted that this question is approached here by treating all 
three variables equally (Feigelson & Babu 1992). It must 
therefore be understood that what is meant in the following 
text by 'departure from the TF relation' is a departure of 
the optimal 'fundamental plane' of spiral galaxy disks from 
the 'TF plane' when all three disk parameters are treated 
equally. Relations constructed for the purpose of galaxy dis- 
tance determinations (see, for e.g., Triay, Lachieze-Rey & 
Rauzy 1994), where, typically, the residuals of only one vari- 
able about some relation are minimised are different; the re- 
sults of such studies cannot be readily compared with the 
current investigation. Consequently, a number of such pub- 
lished data sets are reanalysed here in a consistent fashion. 
The tilt to the optimal 'fundamental plane' of spiral galaxy 
disks is subsequently used to infer the mass-to-light ratio as 
a function of disk scale-length and surface brightness. 

Section 2 introduces the field galaxy samples that have 
been used, and the techniques employed to obtain their disk 
parameters. Section 3 describes the two-dimensional plane 
which is used to represent the location of galaxy disks in the 
three-dimensional space of observables, and describes how to 
relate this plane to the M/L ratio. The results and analysis 
of fitting such a plane to the data sample in hand is presented 
in section 4. The same method of analysis is then applied 
to various literature data samples, and extended to longer 
wavelengths through the study of a homogeneous sample 
of Cluster spiral galaxies for which B, R, I, and K' data is 
available. Lastly, a brief inspection of the M/L ratio depen- 
dency on disk colour is performed. The main results from 
this study are given in Section 6. 



2 PHOTOMETRIC AND SPECTROSCOPIC 
DATA 

Three photometric samples have been compiled, containing 
galaxy profiles in both B- and i?-bands. All these galaxies 
were then structurally modelled using the same algorithm. 
In order to cover a large range of disk parameters, the well- 
defined sample of morphologically undisturbed HSB galaxies 
from de Jong & van der Kruit (1994) and de Jong (1996) 
was used together with two large samples of LSB galaxies 
for which CCD data are available. The first of which is the 
sample of 21 disk-dominated LSB spiral galaxies from de 
Blok et al. (1995). The second is the sample of 21 bulge- 
dominated LSB spiral galaxies from Beijersbergen, de Blok 
& van der Hulst (1999). 

The selection criteria are discussed at length in each pa- 
per, and they have been previously summarised in Graham 
& de Blok (2001). Briefly, de Jong & van der Kruit (1994) 
constructed a statistically complete, diameter- and volume- 
limited sample of 86 HSB disk galaxies having a red minor- 
to-major ratio greater than 0.625 (i.e. inclinations less than 



~50°), and a red major-axis >2.0'. Larger and brighter disk 
galaxies are of course over-represented in such a sample as 
they can be detected at greater distances; equivalently, small 
faint disks will be under-represented. Following McGaugh & 
Blok's (1998) study of LSB galaxies and the TF relation, the 
current objective is simply to expand the baseline of observ- 
able parameters which are fed into such an analysis. This 
HSB sample is therefore supplemented here with two LSB 
galaxy samples; however, as it turns out, the main conclu- 
sions from this study are supported by the analysis of the 
HSB data alone. 

The first LSB sample is from de Blok et al. (1995), con- 
sisting of late-type LSB galaxies from van der Hulst et al. 
(1993) and the lists by Schombert & Bothun (1988) and 
Schombert et al (1992). De Blok et al. randomly selected a 
sub-set of galaxies with inclinations less than 60°, central 
surface brightnesses fainter than 23 B-mag arcsec -2 , and 
having single-dish Hi observations available. The LSB sam- 
ple from Beijersbergen et al. (1999) is a random selection of 
'bulge-dominated' LSB disk galaxies imaged in the ESO-LV 
catalog (Lauberts & Valentijn 1989). Galaxies were selected 
if they had: inclinations less than 50°, diameters of the 26 
B-mag arcsec -2 isophote greater than 1' and smaller than 
3' (chip size limitation) , and a surface brightness of the disk 
at the half-light radii fainter than 23.8 B-mag arcsec -2 . The 
combined sample of HSB and LSB galaxies is not complete 
in any volume-limited sense, but does span more than 4 mag 
arcsec -2 in central disk surface brightness. 

Every profile has been modelled by simultaneously fit- 
ting a seeing-convolved Sersic (1968) r 1//n bulge and a 
seeing-convolved exponential disk. Uncertainty in the rel- 
ative sky-background level can be a significant source of er- 
ror when determining the disk parameters; consequently, the 
parameters for the LSB galaxies (~ 1/3 of the final sample) 
are less secure than those for the HSB galaxies. The LSB 
galaxies have typical errors of ~20 per cent in their disk 
scale-length determinations, while the HSB galaxies gener- 
ally have errors of less than 10 per cent (and typically 7 per 
cent). The profile derived central surface brightness mea- 
surements are accurate to typically 0.1 mag arcsec -2 for all 
galaxy types. They have also all been homogeneously cor- 
rected for: Galactic extinction (Schlegel, Finkbeiner & Davis 
1998), ^-correction (using the tables of Poggianti 1997), 
(1+z) 4 redshift dimming, and inclination corrected follow- 
ing the prescription Tully & Verheijen (1997) derived for the 
Ursa Major galaxy Cluster. 

The inclination of each galaxy (determined from the R- 
band images) have been taken from their respective papers, 
with the exception that the inclinations given in de Blok, 
McGaugh & van der Hulst (1996) were used for the disk- 
dominated LSB galaxy sample. The galaxies have low to 
intermediate inclinations (i<60°). On the one hand, this is 
good because it means that the required surface brightness 
inclination corrections (to a face-on orientation) are small. 
Failing to account for such internal extinction can result in 
substantial errors - these can be greater than 2 mag for 
highly inclined disks. On the other hand, it means that the 
inclination corrections (to an edge-on orientation) for the 
rotational velocity terms will be large. 

The 21 cm velocity data has come from a number of 
sources. The linewidth at 20 per cent of the peak HI inten- 
sity, denoted by W20 and usually measured to an accuracy 
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of a few percent, has been used in all cases. Following Zwaan 
et al. (1995), these values were corrected for both inclination 
and random motion effects, as detailed in Tully & Fouque 
(1985), to give the rotational velocity u rot . Dynamical data 
for the disk-dominated LSB galaxy sample has come from 
de Blok et al. (1996) and van der Hulst et al. (1993), and 
is thoroughly described in these papers. The data for the 
bulge-dominated LSB galaxies has come from de Blok (2001, 
priv. comm.) and was reduced in the same way as the disk- 
dominated sample of LSB galaxies. Measurements of W20 for 
the HSB galaxy sample have come from either the HI cata- 
log of Huchtmeier (2000) or the database of LEDA. Uncer- 
tainties in the true inclination of each disk (due to possible 
intrinsic disk elongation) is likely to be the dominant source 
of error in the rotational velocity measurements and will 
contribute to scatter in the following diagrams. Adding an 
error estimate of 10° to every galaxies photometrically de- 
rived inclination angle resulted in a mean change of 0.09 dex 
in the derived rotational velocities. If more accurate inclina- 
tions derived from kinematical data (see, e.g., Franx & de 
Zeeuw 1992; Andersen et al. 2001) were available, one could 
expect the scatter in the following relations to be reduced. 
Similarly, rotation curve asymmetries have not been dealt 
with (Richter & Sancisi 1994; Haynes et al. 1994; Matthews 
& Gallagher 2002). 

Considering galaxies with both photometric and kine- 
matic data left a sample of 65 HSB and 30 LSB galaxies in 
the B-band, and 64 HSB and 31 LSB galaxies in the _R-band 
(see Table 1). The galaxy F564-v3, the closest galaxy in the 
sample of disk-dominated LSB galaxies - at a distance of 
only 6 Mpc - is the only galaxy excluded at this point due 
to the uncertainty on its true distanced. 



t Inclusion of F564-v3 in the subsequent analysis reveals that 
it lies far from the relation defined by the other 94 galaxies in 
the sample; suggesting its distance is not well indicated by its 
redshift. 
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Table 1. Basic galaxy data 
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Redshifts for eso-lv 035011C and eso-lv 0050050, from 



the bulge-dominated sample of LSB galaxies, were obtained 
from the Parkes telescope HIPASS database (Barnes et al. 
2001). Redshifts for the other galaxies have been taken from 
the NASA/IPAC Extragalactic Database (NED). The he- 
liocentric redshifts were uniformly corrected for Galactic 
rotation, Virgo-infall, infall toward the 'Great Attractor', 
and the 'local anomaly' using the model from Burstein et 
al. (1989). The disk scale-lengths were then converted from 
arcseconds to kpc using a Hubble constant Ho — 75 km s _1 
Mpc -1 . 

In passing, it is noted that the sample of galaxies has 
been taken from the field, rather than from a single cluster, 
and so one can expect to have a greater degree of scatter 
in the relationship between log?;, po, and log/i than might 
otherwise be attained. This is because peculiar velocity flows 
leading to distance errors will introduce errors into the es- 
timates of the scale- lengths. Additionally, like most studies 
to date, W20 has been used to derive v Iot rather than using 
Vmax or Ufl a t which are not available. It is acknowledged that 
the available kinematic data may therefore not be the opti- 
mal estimate to the rotational velocity, although McGaugh 
& de Blok (1998) state that W20 gives a good approximation 

tO tJflat. 



3 THE RELATION BETWEEN logu, ^ AND 
log h, AND IMPLICATIONS FOR THE M/L 
RATIO 

From the balance of forces between centripetal acceleration 
and gravity, one has the dynamical mass M(r) = v 2 r/G, 
where G is the gravitational constant. Combining the previ- 
ous expressions for M and L (section 1) gives 

M v 2 r 

where M and v are measured at the same radius r. A com- 
mon practice is to take r as some multiple of h (e.g. Ah, 
McGaugh & de Blok 1998) and assume the estimate of v 
derived from W20 is a good estimate of not only Ufl at (Mc- 
Gaugh et al. 2000) but also of the rotational velocity at Ah. 
Doing this, one arrives at the virial theorem 



v oc (M/L)I h. 



(2) 



Next, from the measured data of many galaxies, one can 
construct - given it exists - the 'plane' 



v oc Ioh v 



(•3) 



where x and y are the optimal exponents dictating the rela- 
tionship between the disk's structural and dynamical terms. 
Combining equations ^ and [] to eliminate the velocity term, 
one quickly arrives at the relation 



M 

2.5 log — oc (2z 

Lj 



1)/j,o -2.5(22/- 1) log h. 



(4) 



Subsequently, lines of constant M/L will have a slope of 
2.5(2jy — l)/(2x— 1) in the /xo-log h diagram (e.g. Graham & 
de Blok 2001). One can similarly derive that lines of constant 
velocity have a slope of 2.5y/x. 

If the value of v ro t derived from W20, or Ufl at , is an 
appropriate velocity to use at some constant multiple of h 



(e.g. Ah), and galaxy disks are virialised, then the M/L ratio 
derived in equation Q should be proportional to the M/L 
ratio derived from 



M/L | 4h oc v 2 Ah/(I h 2 



(5) 



I shall refer to this as the 'kinematical' or 'dynamical' mass- 
to-light ratio, and that in equation 4 as the 'photometric' 
mass-to-light ratio. 

What M/L ratio does one really want to measure? Lets 
start by looking at the different ways which can be used to 
measure L. The flux within some limiting isophotal radius is 
clearly inappropriate. This is strikingly evident in the case 
of LSB galaxies whose central surface brightness values may 
be fainter than the limiting isophotes used to measure the 
HSB galaxies. In such an instance, the flux measurement of 
the LSB galaxy would be zero. The use of the outer most 
observed radius is also not ideal because it depends on obser- 
vational details such as exposure times, telescope aperture, 
and what signal-to-noise one is willing to accept; that is, the 
outer most observed radius is not an intrinsic quantity of a 
galaxy. Another way to estimate the luminosity of galaxies 
is to measure their flux within some fixed aperture size, say 
20 kpc. Although this is indeed a standard and consistent 
approach, it is again not a favourable one. The reason is 
because, depending on [io and h, this metric will measure a 
different fraction of the total galaxy luminosity. Given that 
exponential disks are structurally homologous, scaling with 
the central intensity and scale-length, one can, and indeed 
should, use these quantities to measure the same fraction of 
disk light for all galaxies (for example, 91 per cent of the 
light from an exponential profile which extends to infinity is 
contained within the inner four scale- lengths; 99.1 per cent 
within the inner 4 effective radii). Therefore, some constant 
multiple of Ioh 2 will represent the same fraction of the total 
extrapolated disk light for all disks. Disk truncation in real 
galaxies at 4 or more scale-lengths will introduce less than 
a 10 per cent error to this approach. 

What about the estimate for the mass Ml One could 
also measure the galaxy mass within 20 kpc, or within a 
fixed number of disk scale- lengths. But is this really what 
we want? If the total (stars plus gas plus dark matter) mass 
profile does not scale with the luminous disk scale- length, 
then estimates for the galaxy mass within a fixed number of 
disk scale-lengths will not correspond to equal fractions of 
the mass distribution. Therefore, a more fundamental quan- 
tity might be the mass within some constant multiple of 
the core-radii of the total mass profile, defined to be the ra- 
dius where the density has dropped by a factor of two from 
its central value. If the mass profiles are homologous, then 
the mass would simply scale with some measure of velocity 
and some measure of radius. The presence of the flat ro- 
tation curves suggests a velocity scale we can use, namely 
Wflat- This then leaves the radial scale-size to be determined, 
and one typically assumes that the scale-size 'a' of the total 
mass profile, scales with the disk scale-length (McGaugh & 

de Blok 1998)0. If they do not, then one should add the term 
— 2.5 log(a/h) to the right hand side of equation 0. 



t Alternative scales could be the velocity vo at the radius ro from 
the density profiles used in Kravtsov et al. (1998). 
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Table 1. - continued 



Galaxy 


Morph. 


Dist 


b/a 


Mo,s 


log(ft s ) 




1 /J, \ 


W 2 o 


Name 


type (T) 


(Mpc) 




— 2 

mag arcsec 


(kpc) 


— 2 

mag arcsec 


(kpc) 


1 —1 

km s 


UGC 12732 


9 


11.4 


0.89 






21.515 


0.087 


127 


UGC 12754 


6 


11.5 


0.70 


21.231 


0.370 


20.207 


0.321 


206 


UGC 12776 


3 


64.1 


0.88 


23.305 


1.355 


21.834 


1.149 


256 


UGC 12808 


3 


54.4 


0.93 


19.839 


0.503 


18.659 


0.504 


362 


F561-1 


9 


70.1 


0.91 


22.738 


0.478 


21.843 


0.437 


069 


F563-1 


9 


53.9 


0.91 


23.716 


0.730 


22.602 


0.628 


115 


F563-V1 


10 


59.5 


0.500 


23.669 


0.450 


22.358 


0.432 


097 


F563-V2 


10 


64.5 


0.88 


22.061 


0.368 






132 


F564-v3 


10 


05.9 


0.87 


23.824 


-0.519 


22.890 


-0.588 


065 


F565-V2 


10 


56.7 


0.50 


24.065 


0.463 


23.394 


0.438 


126 


F567-2 


9 


83.5 


0.94 


23.995 


0.692 


22.936 


0.577 


074 


F568-1 


5 


94.1 


0.90 


23.617 


0.775 


22.619 


0.710 


150 


F568-3 


7 


86.1 


0.77 


22.571 


0.575 


21.611 


0.528 


180 


F568-V1 


7 


84.5 


0.77 


22.741 


0.484 


21.942 


0.466 


176 


F571-5 


9 


65.1 


0.77 


23.362 


0.496 


22.647 


0.448 


179 


F571-V1 


7 


84.2 


0.82 


23.590 


0.511 


22.637 


0.472 


084 


F5741- 


7 


99.2 


0.42 






22.386 


0.749 


223 


F574-2 


9 


91.9 


0.88 


24.102 


0.754 


23.027 


0.635 


068 


F577-V1 


7 


103.7 


0.82 


23.787 


0.799 


23.046 


0.718 


094 


F583-1 


9 


35.1 


0.45 


23.012 


0.318 


22.234 


0.275 


183 


UGC 1230 


9 


50.4 


0.94 


23.208 


0.677 


22.333 


0.609 


115 


UGC 128 


8 


58.9 


0.55 


23.405 


0.864 


22.488 


0.764 


214 


UGC 6614 


1 


92.1 


0.59 






22.567 


1.114 


262 


1150280 


5 


88.2 


0.62 


20.772 


0.657 


19.639 


0.640 


190 


1530170 


1 


88.3 


0.71 






21.369 


0.941 


198 


2520100 


3 


130.9 


0.76 


23.894 


1.406 


21.106 


0.929 


434 


0350110 


1 


86.1 


0.77 


20.829 


0.690 


19.529 


0.648 


285 


3740090 


1 


41.8 


0.81 


22.704 


0.428 


21.145 


0.374 


076 


4220090 


7 


65.9 


0.91 


22.901 


0.717 


21.796 


0.589 


078 


4250180 


7 


95.1 


0.84 


23.378 


0.945 


21.932 


0.821 


101 


4990110 


7 


41.0 


0.83 


22.738 


0.423 


21.995 


0.498 


090 


0050050 


8 


77.2 


0.92 


22.453 


0.705 


21.454 


0.621 


132 


0540240 


2 


198.7 


0.84 


23.555 


1.111 


22.289 


1.012 


399 


5520190 


2 


159.5 


0.62 


22.488 


1.115 


20.810 


0.942 


459 


0590090 


8 


18.4 


0.88 


21.847 


0.283 


20.812 


0.262 


053 


1220040 


1 


127.0 


0.63 


23.272 


1.227 


21.213 


0.882 


101 



Column 1 gives the galaxies identification. The HSB sample from de Jong & van der Kruit (1994) are des- 



ignated according to their UGC number; the LSB galaxies from de Blok et al. (1995) which are labelled 
F (except for 3 galaxies labelled with a UGC number) are from the lists of Schombert et al. (1992); the 
LSB galaxies from Bcijcrsbergen et al. (1999) are denoted simply with a numeric code which refers to the 
ESO-LV number from the catalog of Lauberts & Valentijn (1989). The galaxy type, or numerical stage index 
T (de Vaucouleurs et al. 1991), is given column 2. Column 3 gives the galaxy distance in Mpc (assuming 
Ho = 75 km s _1 Mpc -1 ) after correcting the observed redshift (courtesy of NED) for Galactic rotation, 
Virgo-infall and infall towards the 'Great Attractor' using the model and code from Burstein et al. (1989). 
Column 4 gives the (ii-band) axis-ratio of each disk. The central £?-band disk surface brightness and disk 
scale-length are given in Columns 5 and 6. The surface brightness term has been corrected for the various 
effects described in Section 2 of the text. Column 7 and 8 are the ij-band equivalent of Column 5 and 6. 
Column 9 gives the observed linewidth at 20 per cent of the peak HI intensity. 



Whether or not galaxy rotation curves are homologous 
appears to be something of an unresolved question (van den 
Bosch et al. 2000). Swaters, Madore & Trewhella (2000; their 
figure 4) present H a rotation curves for a small sample of 
LSB and HSB galaxies which appear to overlay each other 
exactly (once normalised with ufl at and h). This contra- 
dicts the work of de Blok & McGaugh (1997) and suggests 
that a/h may be constant. Most recently, fitting pseudo- 
isothermal mass models (in preference to Navarro, Frenk 
& White (1996) CDM halos), de Blok, McGaugh & Rubin 



(2001; their figure 10) studied a sample of well-resolved LSB 
galaxy rotation curves and found no appreciable evidence for 
any correlation between a/h and the maximum rotational 
velocity, or equivalently luminosity. 

In what follows, the additional term a/h will be ignored 
and the kinematical M/L ratio shall be that within 4 scale- 
lengths, as defined in equation [|. 

If the TF relation is indeed the optimal representation 
of disk galaxies, then one must find from equation [] that 
y=2x and v oc L x . The exact slope of the TF relation is im- 
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portant because it tells us the M/L dependency on the other 
galaxy parameters. For example, if L tx v 2 (x=0.5, y—1.0) 
then (M/L) tx h and is independent of surface brightness.pl 
However, if L cx v A (a;=0.25, y=0.5) then (M/L) oc 1/y/h, 
independent of galaxy size. 

In the original paper of Tully & Fisher (1977), they 
found a slope of 2.5±0.3 using photographic magnitudes 
and the velocity line width at 20 per cent of the max- 
imum intensity. More recently, Zwaan et al. (1995) ob- 
tained a B-band slope of 2.6, implying that log(Af/L_g) oc 
(0.09/io+0.54 log h)P\ This result is not at odds with the con- 

' — ' 1/2 

elusion of Zwaan et al. (1995) that M/L oc / , because 
their statement only holds valid, as they state, for galaxies 
of a fixed/constant luminosity. This point has, however, not 
always been appreciated in the literature. 

Let us consider the case where L tx v 2 and log(Af/L) oc 
log h, independent of fio- If one was to sample different galax- 
ies having the same h but different fio, \og(M/L) would show 
no dependency at all on fiQ. If one was to then sample a dif- 
ferent set of galaxies having equal luminosity, but with dif- 
ferent no and h, a plot of log(M/L) versus /to would of course 
show a strong correlation This is because of the relation be- 
tween fio and h at constant L. One must therefore be careful 
of misleading results of this kind. Completely independently 
of any observational results, one can see this degeneracy 
at constant L by dividing the relation M oc v 2 h by L to 
give (M/L) oc v 2 h/L. If the TF relation holds true with any 
slope at all, then one has that (M/L) oc h for any fixed L. Of 
course, for different values of L the constant of proportion- 
ality between M/L and h is different (unless L oc v 2 ) and so 
one cannot use this expression to compare the M/L ratios be- 
tween galaxies of different luminosity, only between galaxies 
of the same luminosity. Alternatively, Zwaan et al. (1995) 
showed that Io (M/L) 2 is constant at any fixed linewidth 
(or, from the TF relation, any fixed luminosity); therefore, 
(M/L) oc 1/vTo- However, if Ioh 2 is constant, one again 
has that (M/L) oc h. Identically, one cannot use the rela- 
tion (M/L) oc l/y/To to compare the M/L ratio of galaxies 
at different luminosities (unless L oc v 4 ). Only for galaxies 
having the same fixed luminosity can these expressions be 
used to indicate differences in the M/L ratio. How the M/L 
ratio of galaxies with different luminosities varies (simulta- 
neously) with scale-length and the central disk intensity is 
what this study shall investigate. 



4 DATA ANALYSIS AND RESULTS 

A Principle Component Analysis was performed on the data, 
revealing that 95 per cent (B-band) and 94 per cent (R- 
band) of the variance within the 3-dimensional space of 
variables (fio, log h, log?;) resides within a two-dimensional 
plane[Jj From this it is concluded that it is indeed appro- 
priate to represent the data with a two-dimensional plane. 



S This result actually holds for any value of y as long as 2=0.5. 
" n.b. Along lines of constant luminosity in the /to - log h diagram, 
if log(M/L) oc 0.1/to + 0.5 log h then it implies that both /to and 
logh contribute equally to changes in M/L. 

II The principle eigenvector accounted for 49 per cent (B-band) 
and 55 per cent (R-band) of the variance. 
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b=-0.36 
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Figure 1. Correlation between the dynamical and structural 
properties of spiral galaxy disks. The 'mixing' value b (see, e.g., 
Djorgovski & Davis 1987) is chosen such that it yields the highest 
linear correlation coefficient between log v Ia t and (log h + 6/io ) • 
The solid line represents the best fit from an orthogonal regression 
analysis, and the dashed lines mark the boundary which is a factor 
of 2 away in v TO t from the line of regression. The different symbols 
denote from which galaxy sample the data are from: filled circles 
represent galaxies from the sample of de Jong & van der Kruit 
(1994), stars are for galaxies from de Blok et al. (1995), open 
circles are for galaxies from Beijersbergen et al. (1999), and the 
plus signs denote the five outlying points which have been rejected 
from the analysis. The perpendicular scatter in both panels a) B- 
band and b) i?-band is 0.11 dex. 



In so doing, five galaxies (3 LSB and 2 HSB) from the fi- 
nal sample of 94 ha ve been identifi ed as potential outliers 
- one in particular (eso-lv 1220040). The following analy- 
sis has therefore been performed excluding |eso-lv 122004C , 
and then performed again excluding all five galaxies, which, 
in a number of tests, consistently lied away from the main 
data. These points are shown in Figure [j] - identified as ly- 
ing more than a factor of 2 in v ro t from the best-fitting line 
to the data sample. Twop^| of these galaxies have an incli- 
nation of ~25°, and another is the closest galaxy from the 



HSB sample - having a redshift of only 800 km s~ . The 
above percentages increased by 1 per cent, to 96 per cent 
(B-band) and 95 per cent (B-band), when all five outlying 
data points were excluded. 

To construct the best-fitting plane, a three-dimensional 
bisector method of linear-regression was used (Graham & 
Colless 1997). The advantages and applicability in treating 
all the variables equally in this manner, as opposed to a 
simple linear regression of one variable on the others, is elu- 
cidated in Isobe et al. (1990) and Feigelson & Babu (1992). 
Basically, one wants to treat all the data/variables in an 
equivalent manner when performing a comparison with the- 
ory. A simple ordinary linear regression (OLS) of one vari- 
able on the other two variables would not achieve this goal, 
and would produce different results depending on which vari- 
able one regressed upon the others. Of those methods which 
do provide a symmetrical treatment of the data, the bisec- 
tor method provides the smallest variance, or uncertainty, 
to the estimated values of the fitted plane. 

Performing this regression, the exponents x and y 



** An additional 10 galaxies have inclinations <25°, but all 
closely follow the relation denned by the other galaxies. 
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Figure 1 — continued 
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Figure 2. Same as Figure [l]a and ^3, except that the value of 
b has been fixed equal to -0.2 to mimic the TF relation, where 
L oc v a and a is obtained here by dividing 2.0 by the slope of the 
fit. The perpendicular scatter in both panels a) B-band and b) 
i?-band is 0.11 dex. 

from equation ^| were computed to be 0.53±0.05 and 
0.79±0.09 (B-band) and 0.45±0.04 and 0.71±0.09 (R- 
band). Therefore, (M/L B ) oc jo.osio.io^o.ssio.is ( B _ band ) 



and (M/Lr) oc 7, 



-0.10±0.08 r0.42±0.18 



(i?-band). Lines of con- 



stant velocity subsequently have a slope of 2.5y/x = 3.7±0.5 
(B-band) and 3.9±0.6 (i?-band) in the fio~\og h diagram. 
A plot of logv against log(lQh y ) does, as one would hope, 
have a slope of 1. This figure is not presented here be- 
cause it is attained by simply dividing the :r-axis in Fig- 
ure |l| by the slope of that data. Therefore, a scaled version 
of such a plot is already displayed. Both of the TF rela- 
tions shown in Figure H (obtained by setting the 'mixing' 
value b = —x/(2.5y) = —0.2) have a vertical scatter of 0.14 
dex, slightly greater than the value 0.13 dex (c.f. 0.65 mag, 
Kannappan, Fabricant & Franx 2002) obtained in Figure hi 
The above analysis was repeated, removing those 
additional four data points marked as outliers in Fig- 
ure [j] This time z=0.50±0.05 and y=0.77±0.07 (B- 
band), and x=0.43±0.04 and y=0.69±0.07 (ii-band). Here, 
(M/L B ) oc /0-°°±°-i° ft o.54±o.i4 (B_band) and (M/L R ) oc 
jr-o.i4iQ.o8jjO.38io.i4 (#_band), and lines of constant veloc- 
ity have a slope of 3.9±0.5 (B-band) and 4.0±0.6 (i?-band) 
in the /io~log h diagram. At least in the B-band, the data 
reveal that the M/L ratio is in fact largely independent of 
surface brightness, and roughly varies as yfh. Lower surface 



0.4 




b/a 



Figure 3. The vertical residuals about the best-fitting line in 
Figure [l]a are plotted against a) each galaxies observed minor-to- 
major axis ratio (b/a), and b) the B — R colour of each disk. 



brightness does not necessarily imply a higher M/L ratio. 
The residuals are shown not to correlate with either (B — R) 
colour (Figure []), or galaxy ellipticity, and hence no obvi- 
ous additional dependency on inclination is expected. Stud- 
ies which explore the density profiles of dark matter halos 
associated with LSB galaxies should perhaps reconsider the 
assumption that all of these systems are dark matter domi- 
nated. 
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Figure 4. The logarithm of the difference between the mass- 
to-light ratio estimates obtained using equation |lj and equation ^ 
have been plotted against a) the central disk surface brightness, 
b) the logarithm of the disk scale-length, and c) the absolute 
magnitude of the disk. No strong trends exists. These plots show 
the B-band data, excluding the 5 outlying points identified in 
Figure |l| The values for x and y used in equation ^ have been 
taken to be 0.50 and 0.77. Constants of proportionality have been 
ignored; the constant c in the x-axis is the same for all galaxies. 



Differences in the mass-to-light ratios estimated from 
equation ^ and equation ^ may signify that galaxies are ei- 
ther not virialised or that the mass profiles do not scale sim- 
ply with the luminous exponential disk scale-length (section 
3). The logarithm of the ratio of the 'kinematic' and 'photo- 
metric' mass-to-light ratios has been plotted against the cen- 
tral disk surface brightness, the logarithm of the luminous 
disk scale- length, and the logarithm of the disk luminosity in 
Figure ^[ No strong correlations exist, although a weak cor- 
relation (Pearson's r = 0.4) with absolute magnitude does 
exist at the 99.97 per cent confidence level. This may be in- 
dicating that the ratio of the mass profile scale-size (a) to 
the (B-band) luminous disk scale-length (h) may correlate 
weakly with luminosity, in the sense that a/h decreases with 
increasing disk luminosity. 

One way to test if the above trend is true would be 
from an analysis of actual rotation curves. If rotation curves 
are homologous, scaling in amplitude with Wg at and radially 
with disk scale-length h, then they should all overlap each 
other once normalised in this way. However, because rota- 
tion curves can rise slowly as they approach Ufl at (see, for 
e.g., Rubin et al. 1985), it is difficult to quantify at exactly 
which radius Ufl at is reached. It can therefore be difficult to 
gauge if rotation curves are homologous. It is proposed here 
that one might be able to use the radius where the velocity 
equals one-half of Ufl at as a measure of the radial scale for 
the velocity profile. Due to the relative steepness of rotation 
curves around this value (after they have been scaled with 
the disk scale- length), errors in vnat should have a reduced 
effect on estimates of this radius. A plot of this radius, in 
units of disk scale-length, against disk luminosity may reveal 
if the assumption that the mass profiles scale with the lumi- 
nous disk profiles is valid, or if the slight correlation found 
above is real. 



4.1 Further analysis 

A range of methods have been used in the literature for 
the task of constructing a two-dimensional plane to three- 
dimensional data. Given the lack of universal agreement as 
to which is the best approach to take, a second (different) 
approach is employed here. Following Djorgovski & Davis 
(1987), the 'mixing' value b that gave the highest linear cor- 
relation coefficient between log v and (log h + 6/io) was com- 
puted. A line of regression to this two-parameter data set 
(once the value of b had been determined) was then fitted; 
in fact, six different lines of regression were fitted using the 
program SLOPES from Feigelson & Babu (1992). The re- 
sults are given in Table ^, and were obtained excluding the 
1 outlying data point ( aso-lv 1220040| ). Through a bootstrap 
resampling of the data an estimate for the error on b was 
obtained and is also given in Table ^| When all five outlying 
data points were excluded, &=— 0.256 ± 0.030 (B-band) and 
b—— 0.265 ± 0.036 (i?-band), in agreement with the results 
obtained in Table bl Again, the analysis and data suggests 
that in the _B-band x ~ 0.5 and the M/Lb ratio is largely 
independent of surface brightness. Once more, only mild in- 
consistencies with the TF expectation value of b — —0.2 
were found using this alternative technique to compute the 
best-fitting plane. Within the 2a confidence interval, the op- 
timal values of b are consistent with the value -0.20. Forcing 
the value of b to equal -0.20 introduced a slight dependency 
of M/L on the surface brightness term - as evidenced by 
the departure of x from the value 0.5 in Table ^. The extent 
of this effect was greater in the redder passband. 

A significant fraction (1/2) of the LSB galaxies have ro- 
tational velocity estimates less than 100 km s _1 . In order to 
obtain an estimate as to how the LSB galaxies may be in- 
fluencing the above results, the analysis was repeated using 
the HSB galaxy sample on their own, and then again using 
only the LSB galaxy sample. Excluding the outliers identi- 
fied in Figure [j], the three-dimensional bisector method of 
regression for the HSB galaxy sample [N = 62, 61) gave v tx 

/ 0.60±0.12 /l 0.82±0.12 (jB _ band) and „ x j0.50i0.09 ft 0.66±0.09 

(J?-band)fT] 

One can therefore conclude that the LSB galax- 
ies are not responsible for the observed, albeit slight, depar- 
tures from the TF relation, nor are they biasing the analysis 
to produce values of x as high as 0.5 and thereby caus- 
ing the M/L ratio to appear largely independent of sur- 
face brightness. The regression analysis of the LSB galaxy 
sample (N = 25,27) gave v oc /°-58±o.i7 ft o.92±o.r 6 ( B _ band ) 
and v oc / -,6±o.i7 /l o.79±o.i8 ( R _ hand ^ fmly consistent with 

the result from the HSB galaxy sample. The reason for the 
slightly larger exponents when the HSB and the LSB galax- 
ies are considered on their own, as opposed to taken to- 
gether, can be easily understood by looking at the distribu- 
tion of LSB and HSB galaxies in Figure hi This slight dif- 
ference in slopes also illustrates how a selection bias against 
LSB galaxies would operate; removing the left half of the 
data in Figure |l| results in an increased slope in the best- 
fit to the remaining data (see, e.g., Lynden-Bell et al. 1988, 
Appendix B). 



tt From the second method of analysis, the 'mixing' values were 
determined to be 6=0.244±0.044 (B-band) and fe=0.260±0.050 
(i?-band). 
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Table 2. Regression analysis: Field Galaxies 

method intercept slope (= y) 



-2.5 X b X y 



B-band: 



-0.265 ± 0.034 



OLS(Y|X) 
OLS(X|Y) 
OLS bisector 
Orthogonal 
Reduced major axis 
Mean OLS 



5.149±0.267 
7.158±0.417 
6.060±0.270 
5.719±0.330 
6.024±0.276 
6.153±0.286 



0.576±0.052 
0.967±0.081 
0.753±0.052 
0.687±0.064 
0.747±0.054 
0.772±0.056 



0.382±0.053 
0.641±0.086 
0.499±0.063 
0.455±0.064 
0.495±0.063 
0.511±0.065 



B-band: fe=-0.200 

OLS(Y|X) 

OLS(X|Y) 

OLS bisector 

Orthogonal 

Reduced major axis 

Mean OLS 



4.539±0.240 
6.321±0.347 
5.326±0.225 
5.134±0.298 
5.305±0.234 
5.430±0.238 



0.633±0.065 
1.112±0.093 
0.844±0.060 
0.793±0.080 
0.839±0.063 
0.872±0.064 



0.317±0.033 
0.556±0.047 
0.422±0.030 
0.397±0.040 
0.420±0.032 
0.436±0.032 



R-band: 6=-0.278 ± 0.040 

OLS(Y|X) 4.754±0.235 0.499±0.046 0.347±0.055 

OLS(X|Y) 6.434±0.360 0.825±0.071 0.573±0.089 

OLS bisector 5.532±0.244 0.650±0.048 0.452±0.067 

Orthogonal 5.127±0.284 0.571±0.056 0.397±0.065 

Reduced major axis 5.489±0.248 0.642±0.049 0.446±0.067 

Mean OLS 5.594±0.256 0.662±0.050 0.460±0.069 



i?-band: 6=-0.200 

OLS(Y|X) 

OLS(X|Y) 

OLS bisector 

Orthogonal 

Reduced major axis 

Mean OLS 



4.230±0.208 
5.727±0.294 
4.902±0.201 
4.663±0.258 
4.876±0.208 
4.978±0.211 



0.580±0.059 
1.003±0.084 
0.770±0.058 
0.702±0.074 
0.762±0.060 
0.791±0.060 



0.290±0.030 
0.502±0.042 
0.385±0.029 
0.351±0.037 
0.381±0.030 
0.396±0.030 



Regression analysis betw een log v and (l og h + bfio) for the field galaxy data presented 



in Section 2 (excludin. 
plane such that v oc I£ h y where fio 
b = -0.2. 



so-lv 122004C). The values x and y define the best-fitting 
2.5 log Iq. The TF relation is obtained when 



A remaining point of concern is that there may be rela- 
tively important distance errors for some of the more nearby 
galaxies due to streaming motions above the smooth Hub- 
ble flow. To address this concern, the regression analysis 
has been repeated excluding the 34 galaxiesp] within 3,000 
km s . The results obtained from the regression analysis 
of this data are fully consistent with the results from the 
larger galaxy sample, with v oc /°-53±o.06 fe o.85±o.i2 ( Bband ) 
and v oc j°-45±o.o5 ft o.75±o.i3 (^.hand). Additionally, correc- 
tion for Virgo-infall alone, using the 220-model of Kraan- 
Korteweg (1986), gave entirely consistent results with those 
obtained using the more detailed velocity flow model of 
Burstein et al. (1989). 

A meaningful comparison of the best-fitting planes for 
the early-type and late-type (>Sc) spiral galaxies in the 
present data set is not possible because the uncertainty on 
the exponents is in general too great. What we can, however, 
conclude at this point is that in the B-band, the near inde- 
pendency of M/Lb on surface brightness holds when the 
HSB galaxies are taken on their own, when the HSB and 
LSB galaxies are taken together, and when one constrains 

M Of the 5 galaxies identified previously as 'outliers' in Figure [l| 
only two of these are within this volume. 



the 'mixing' value b to -0.2 in order to mimic the TF rela- 
tion. Although, it is noted that this situation is expected to 
change in the redder passbands, where the slope of the TF 
relation is steeper (see section 5). 



4.2 Does M/Lb correlate with ^to,s? 

Two points can be made regarding previous B-band studies 
which have claimed, and even plotted, a trend between mass- 
to-light ratio and central disk surface brightness. The first 
point is addressed by plotting the present data. Figure 
shows the 'kinematical' M/Lb ratio (equation^) versus the 
central B-band disk surface brightness [io,B- One can see 
that although an upper envelope appears to exist, the cur- 
rent data present no strong correlation between M/Lb and 
Ho,b- Furthermore, the upper envelope is merely a reflec- 
tion of the upper (bright) envelope in the fj,o~\og h diagram 
(Graham 2001b; Graham & de Blok 2001) which is due to 
the absence of large scale-length galaxies with bright disks. 
Panel b) in Figure ^ displays the data from Table 1 of Mc- 
Gaugh & de Blok (1998). Their tabulated circular velocities, 
B-band scale-lengths and central disk surface brightnesses 
were used to compute the kinematical mass-to-light ratio 
within four scalelengths. This was done in exactly the same 
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Figure 5. The kinematic M/L ratio (in solar units) is shown 
against the central disk surface brightness (B-band data). Panel 
a) contains the data described in Section 2, excluding the five 
outlying points identified in Figure |l|. Panel b) shows the smaller 
data sample from Table 1 of McGaugh & de Blok (1998). Panel 
c) shows the data from panel a), excluding LSB galaxies with /Uo 
fainter than 22.65_B-mag arcsec -2 and with h smaller than 8.45 
kpc (see Section for details). Panel d) shows the HSB galaxy 
sample from panel a) (excluding all three HSB galaxies with fio 
fainter than 23.0 B-mag arcsec~ 2 ) together with the LSB galaxies 
in common with McGaugh & de Blok (1998). 



way as for panel a) and is in fact the method claimed to 
have been used in McGaugh & de Blok (1998). Strangely, 
their tabulated M/L ratios are different to the values de- 
rived from these basic parameters, and consequently their 
Figure 4b looks different to Figure ^(b shown here. Although 
the upper envelope is common between panels a) and b), 
the lower-left portions differ significantly. 

The disk-dominated LSB galaxies from de Blok et al. 
(1995) have roughly the same scale-lengths as normal HSB 
spiral galaxies - their median B-band scale-length is 4.3 kpc 
(Ho = 75 km s -1 Mpc" 1 ; de Blok et al. 1995). The bulge- 
dominated LSB galaxies from Beijersbergen et al. (1999) 
have a larger median Z?-band disk scale-length of 12.6 kpc 
(Beijersbergen et al. 1999). Taking 8.45 kpc (=[12.6+4.3]/2) 
to represent some kind of division between the large and 
small (or rather, large and normal) sized LSB galaxies, 
panel c) shows the results after selecting against (i.e. ex- 
cluding) those disks with central surface brightnesses more 
than one mag fainter than the canonical Freeman (1970) 
value of 21.65 B-mag and with scale-lengths less than 8.45 
kpc. These galaxies should roughly have, on average, the 
same scale-lengths as the HSB galaxy sample, but fainter 
surface brightnesses. If the M/L ratio is independent of 
surface brightness, then we should find that we have re- 
moved those LSB galaxies with similar M/L ratios as the 
HSB galaxies - which indeed is what we do find. Panel d) 
shows the current HSB galaxy sample (excluding the three 
disks with fj,o fainter than 23.0 B-mag arcsec -2 ) together 
with those LSB galaxies in common with McGaugh & de 
Blok (1998). From this restricted sample, a somewhat sim- 
ilar correlation between M/ Lb and fio,B to that present in 
the data of McGaugh & de Blok (1998) appears. The de- 
pendence of the B-band mass-to-light ratio on the central 



disk surface brightness can therefore arise from ones sample 
selection, even with a range of disk luminosities 

McGaugh & de Blok (1998) excluded a number of LSB 
galaxies for a variety of different reasons. All of the LSB 
galaxies in common with both studies but excluded by Mc- 
Gaugh & de Blok (1998) (N = 7) reside very close to the 
mean relation defined by the other galaxies in Figure That 
is, they do not appear to be outliers. 

The second point which should however be mentioned 
is that in contradiction with most B-band TF studies, the 
small sample of McGaugh & de Blok (1998) apparently has 
a TF slope of almost 4, i.e. x = 0.25 (their Figure 1). In 
this case, as explained in Section 3, one does expect to find 
a strong correlation between M/L and /io, as these authors 
did for their galaxy sample. 



5 LITERATURE COMPARISONS 

A somewhat similar three-dimensional analysis of 511 HSB 
Spiral (and Irregular) galaxies having types Sa through Im 
has been performed by Burstein et al. (1997) in the B- 
band. While they also used rotational velocity estimates 
from W20 measurements, the nature of their photometric 
data is different to that used here. Using the RC3 (de Vau- 
couleurs et al. 1991) tables of B-band circular effective aper- 
ture size and mean enclosed surface brightness, Burstein 
et al. (1997) used the axial-ratio-based statistical relation- 
ships in the RC3 to derive the mean effective surface bright- 
ness <fj.> e = — 2.51og</ e > within the effective radius r e . 
They then corrected <fj,> e for Galactic extinction follow- 
ing Burstein & Heiles (1982) and for internal extinction us- 
ing the correction 1.51og(b/a). Recognizing errors of more 
than 0.1 dex in their photometric parameters, Burstein et al. 
(1997) chose not to apply the small (1 + z) 4 correction or any 
K- correction term. Their investigation effectively used the 
total (bulge + disk) galaxy half-light radius and mean sur- 
face brightness within this radius. Depending on the contri- 
bution of light from the bulge, varying fractions of disk scale- 
length, and hence the disks themselves, were consequently 
sampled in their analysis. The present investigation differs 
as it has used one complete disk scale-length and the cen- 
tral disk surface brightness; an approach, it should be noted, 
which effectively regards the bulge matter as dark matter - 
the luminous flux from the bulge is ignored, while the pres- 
ence of the bulge permeates into the analysis through the 
rotational velocity of the disk. 

If the relative contribution from the bulge light^ is 
small enough that it can be neglected and the galaxy treated 
as simply an exponential disk, then the use of r e and </i> e 
(or /He, the surface brightness at r e ) should yield similar re- 
sults as the use of h and no- This is because r e — 1.67/i 
and <n> e = fJ-o + 1-12 (fj, e — (io + 1.82) for an expo- 
nential disk. To investigate the influence of the bulge, and 
perform a somewhat more direct comparison with Burstein 
et al. (1997), the half-light radius and the surface bright- 
ness at this radius have been derived for the present galaxy 

S§ B/D ratios for much of the present galaxy sample have al- 
ready been presented in Graham (2001a, b) and Graham & de 
Blok (2001). 
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Figure 6. a) The B-band half-light radii of the disks (r c disk = 
1.67A) are plotted against the B-band half light radii of the galax- 
ies (f e ,galaxy)- The presence of the bulge accounting for the dif- 
ferences. The dashed line denotes a 20 per cent difference from 
agreement (the solid line). The marked outliers have very promi- 
nent bulges within the half-light radius of the disk, b) The disk 
surface brightness at r disk (= fio + 1.82) is plotted against the 
galaxy surface brightness at r ega i axy . 
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sample and the fundamental plane recomputed using these 
values. The differences between the half-light radius of the 
disk (=1.67ft) and the half-light radius of the galaxy (r e ), 
and the surface brightness of the disk at 1.67/i and of the 
galaxy at r e are shown in Figure |^. With the exception 
of a few galaxie sp^] having very prominent bulges within 
their disk half-light radius, the disk and galaxy half-light 
radius usually varied by less than 20 per cent, and the sur- 
face brightness terms by less than ~0.4 mag arcsec -2 . The 
optimal fundamental plane constructed using r e and /i e for 
the full sample (excluding the previously mentioned outliers) 
is v oc /0-53±o.o5 r o.79±o.o9 (jB _ barid ). For the jjSB sample 
one finds v oc /0.6i±o.u r o.80±o.ii j and for the LSB sample 

v (X J 0.57±0.13 T ,0.87±0.17_ 

From the sample of 386 Sa-Sc Spiral galaxies analysed 
by Burstein et al. (1997), they found v oc< I >0- 57 ± 03 
r o.85±o.05 (73_b and ) Applying the three-dimensional bisec- 
tor regression analysis used in this paper to this very 
same data sample from Burstein et al. (1997) gives v oc< 



^ These galaxies are not the outliers identified in Figure [l]. 



/ >^ 5±0 - M r .o 78±o.os ( B _ band )_ From equation | this im- 
plies that (M/L B ) oc< / >o io ± 08 r o.58±o.i 0j in good agree- 
ment with the present papers investigation. 

When Burstein et al. (1997) included all Spiral and 
Irregular galaxy types in their analysis, they found v oc< 
j > o.39±o.02 r o.6i±o.02 ( B _ band ) ] a result with notab l y lower 

exponents. Although, applying the three-dimensional bisec- 
tor regression analysis gives v oc< I >o 41 ± 01 r o.62±o.02^ 
and excluding the Irregular type galaxies raises the expo- 



nents slightly to give v oc< / > 



0.44±0.02 „0.66±0.02 



(B-band). 



When the same internal extinction correction as used by 
Burstein et al. (1997) is applied to this papers field galaxy 
sample, one obtains v oc /0.48±o.o4 r o.73±o.o7 ( B _ band ) ] con _ 

sistent at the ler level with the result from Burstein et al. 
(1997). 

The result from Burstein et al. (1997), using all 511 disk 
galaxies and irregular galaxies, leads to (although not given 
in their paper) L oc l ,2.56±o.i3 ? .o.44±o.07 ( B _ band ). An even 
larger departure from the TF relation and dependence on 
a third parameter was found by Willick (1999) who derived 
L oc v 2A r ( l' 7 (ii-band). Others have also previously found 
evidence for the residuals about the TF relation to correlate 
with a third parameter (Kodaira 1989; Feast 1994; Koda 
et al. 2000a; Han et al. 2001, Kannappan et al. 2002), and 
it may be that the TF relation reflects closely, but not ex- 
actly, the Fundamental-Plane-like distributions shown here. 
However, there are numerous papers (e.g. Strauss & Willick 
1995; Courteau 1997; and references therein) which, like this 
one, do not detect deviations above the l-2cr level. 

It should also be kept in mind that different papers have 
used different methods of regression analysis to fit a plane or 
line to their data points. The present papers interest lies in 
deriving M/L ratio estimates which can be compared with 
theory. The disk parameters have consequently been treated 
equally when constructing the optimal plane; an approach 
which differs from the numerous direct and inverse TF rela- 
tions constructed in the past for the determination of galaxy 
distances - which is not the subject of the current paper. A 
detailed comparison with every previous study is therefore 
not straight-forward if at all even possible. However, a de- 
tailed analysis of one other well-respected data sample shall 
be presented. 

The reduced and tabulated photometric data of Tully 
et al. (1996) and the reduced kinematic data from Verheijen 
& Sancisi (2001) is analysed. This consists of B, R, I and K' 
scale-lengths and central disk surface brightnesses, obtained 
from marking the disk by eye, and rotational velocities taken 
from where the rotation curve is flat (tifiat)- The galaxy sam- 
ple are all from the Ursa Major Cluster and consequently 
distance errors will not contribute quite as much to the scat- 
ter. For the sake of consistency, the surface brightness term 
has been uniformly corrected here in the same way as the 
previous data was dealt with (see Section 2). Because this 
new sample contains relatively edge-on Spiral galaxies, the 
iT'-band inclination corrections to the central surface bright- 
ness terms are large, in some cases >2.0 mag arcsec -2 . For 
four galaxies which were catalogued to be completely edge- 
on, their inclinations were set to 85° when calculating this 
correction. 

The results of the three-dimensional bisector regres- 
sion analysis are presented in Table |S| In the B-band, the 
analysis reveals that log(Ai/Ls) oc (0.08/io +0.46 log ft), 
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Table 3. Regression Analysis: Cluster Galaxies 
Band x y M/L 



B 
R 
I 

K' 



0.396±0.059 
0.334±0.035 
0.303±0.029 
0.258±0.020 



0.732±0.118 
0.642±0.098 
0.601±0.092 
0.569±0.073 



j-0.208±0.118 ft 0.464±0.236 
j-0.332±0.070 ft 0.284±0.196 
j-0.394±0.058^0.202±0.184 
j-0.484±0.040 ?l 0.138±0.146 

Three-dimensional bisector regression analysis between logti, /ig, 
and log h for Ursa Major Galaxy Cluster members. The exponents 
x and y define the best-fitting plane such that v oc I^h v where 
^0 = — 2.5 log /o- The photometric data (fio,h) comes from Tully 
et al. (1996) and the kinematic data (v) from Verheijen & Sancisi 
(2001). 




K- band (C=1.0) 

slope=0.51 

b=-0.19 



-3.8 -3.6 -3.4 -3.2 -3 -3.8 -2.6 

log(h) + b/x 

Figure 7. Orthogonal regression analysis to the A"'-band data 
from Tully et al. (1996) and the velocity data from Verheijen & 
Sancisi (2001). The value of b was computed, as in Figure pi such 
that it resulted in the highest linear correlation coefficient be- 
tween log2tifl at and (log h + bfio). When b = 0.20 these diagrams 
represent the TF relation, where L oc ^ 2 / slc >P c . The value of C 
dictates the apparent opacity of spiral galaxy disks for determi- 
nations of /iq, such that the inclination correction to the surface 
brightness term fio can be written as 2.5C log(a/6), where a/b 
is the observed outer disk axis ratio. Panel a) assumes the disks 
are opaque (C = 0), while panel b) assumes that they are fully 
transparent (C = 1), resulting in considerably less scatter to the 
diagram. 



in good agreement with the previous B-band results and 
revealing that changes in scale-length are more important 
than changes in surface brightness. Figure m shows the re- 
sults for the K'-band data with (C = 1.0) and without 
(C — 0.0) the inclination correction (2.5Clog[a/6], where 
a/b is the observed disk axis-ratio) to the surface brightness 
term. From the scatter one can clearly see that at near- 
infrared wavelengths, spiral galaxy disks should be taken to 
be transparent for such tests (see also Graham 2001b). One 
can see that the TF relation is obtained (i.e. y ~ 2x), and 
the slope increases from the B-band to the if'-band such 
that Lb oc u 2 ' 5 and L K i oc v 4 ' . The reduction in scatter 
with increasing wavelength has also been observed for many 
years (Aaronson, Huchra & Mould 1979). Therefore, in the 
near-infrared, where one is affected little by obscuring dust 
and young stellar populations, the M/L K i ratio varies as 
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Figure 8. The logarithm of the difference between the mass- 
to-light ratio estimates obtained using equation [ij and equation ^ 
have been plotted against the central disk surface brightness (top 
panel), the logarithm of the disk scale- length (middle panel), and 
the absolute magnitude of the disk (bottom panel). These plots 
show the B-band data from Tully et al. (1996) and the velocity 
data from Verheijen & Sancisi (2001). The values for x and y used 
in equation Q have been taken to be 0.396 and 0.732 (see Tablc|i|). 
Constants of proportionality have been ignored, and the value of 
c in the x-axis is simply a constant. 



To again investigate possible differences between the 
mass-to-light ratios obtained using equation ^ and equa- 
tion and hence possible signatures that the total-mass 
profile does not scale linearly with the luminous disk profile, 
we have repeated Figure |i| this time using the Ursa Major 
galaxy Cluster data. Figure ^ presents the B-band analysis, 
which can be seen to display an identical behavior to that 
observed in Figure This tends to suggest that the ratio 
of the mass profile scale-size to the B-band luminosity disk 
scale-length may not be constant, but scales weakly with 
the B-band disk luminosity (see Section 3) . Figure M shows 
the Jf'-band data where one can see that the weak B-band 
dependency on absolute magnitude (r = 0.4, Figure |j| lower 
panel) is reduced even further. 

Lastly, it should perhaps also be noted that there is an- 
other mechanism to change the slope of the optimal plane 
within a given passband. Tully et al. (1998) have presented 
evidence that the appropriate extinction correction for disk 
galaxies in Clusters may be a function of galaxy luminos- 
ity. This was a subsequent paper to the Tully & Verheijen 
(1997) paper which addressed the issue of extinction correc- 
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Figure 9. Same as Figure |j but using J<"'-band data. The values 
for x and y used in equation Q to derive (M/L) p i lot have been 
taken to be 0.258 and 0.569. 



tions and provided the somewhat more standard corrections 
which have been used in this paper. Tully et al. (1998) anal- 
ysed the data from two galaxy Clusters (53 galaxies from 
Ursa Major and 34 galaxies from Pisces) to determine a 
TF slope of 3.1 in the B-band after applying their luminos- 
ity dependent extinction correction. This is an increase of 
0.2 from 2.9 found by Sprayberry et al.'s (1995) study of 
the Ursa Major Cluster. Perhaps contrary to this, Willick 
et al. (1996) found no evidence for a luminosity-dependent 
extinction correction when using field galaxies. The inter- 
esting analysis by Tully et al. (1998) raises many questions. 
For instance, the data for the Pisces Cluster galaxies on their 
own appear to suggest no luminosity dependent extinction 
correction. Additionally, whether the results in Tully et al. 
(1998) can be extended to LSB galaxies is unclear; why an 
unevolved LSB galaxy with the same luminosity as a HSB 
galaxy should have the same amount of extinction due to 
dust is not obvious. Why the SO and SO/a galaxies do not 
conform is a further mystery. The issue of corrections due 
to dust in field galaxies is not yet fully understood, or at 
least commonly accepted, and therefore it should be noted 
that the extinction corrections used here may need adjusting 
when this subject is on a firmer footing. 
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Figure 10. The logarithm of the dynamical mass-to-light ratio 
obtained using equation 1 (with r = Ah) is plotted against two 
measures of the disk colour (B — R). The left hand column shows 
the colour difference between the integrated disk magnitude in B 
and R, while the right hand column shows the difference between 
the central disk surface brightness in B and R. Only the disk 
luminosity (not the bulge and disk luminosity) has been used to 
derive the mass-to-light ratio, which is given here in solar units. 
An orthogonal regression analysis has been used to derive the 
slope (and error) to the correlation, and the Pearson correlation 
coefficient is also given for those plots displaying a possible cor- 
relation. 



6 COLOURS 

The (B — R) colour of each galaxy's disk has been estimated 
using two approaches. Having previously modelled both the 
bulge and disk, we are better able to determine the (B — R) 
colour for the disk alone - rather than the integrated colour 
of the bulge and disk. The first method has taken the dif- 
ference between the corrected (see section 2) central disk 
surface brightness in B and R, and the second approach 
has used the difference between the total disk magnitude 
determined in the B and R passbands. Each galaxies dy- 
namical mass-to-light ratio estimate (from equation |E|) is 
plotted against its B — R colour in Figure |Io| for the HSB 
and LSB field galaxies, and in Figure hll for the Ursa Major 
Cluster data from Tully et al. (1996) and Verheijen & San- 
cisi (2001). Only the disk light (that is, not the bulge light) 
has been used in this derivation of the dynamical mass-to- 
light ratio. Using the combined disk and bulge light to de- 
termine the mass-to-light ratio has very little effect on this 
result. It is immediately obvious that, within the noise of 
the present field galaxy data, no obvious correlation exists 
between the B- and i?-band dynamical mass-to-light ratio 
and B — R colour. An orthogonal regression analysis has 
been performed for those plots showing a possible correla- 
tion, and the Pearson correlation coefficient determined (re- 
sults shown within the figures). In general, the correlations 
are weak. 

Using a suite of simplified galaxy disk evolution mod- 
els, Bell & de Jong (2001) have presented the case that 
the stellar mass-to-light ratio of the disk (M*,disk/£disk) 
should be correlated with colour. Their strongest correlation 
was obtained with the colour (B — R) such that the loga- 
rithm of a disks B-band stellar mass-to-light ratio varies as 
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Figure 11. The same notation applies here as in Figure [uj. This 
time the Ursa Major galaxy Cluster data introduced in section 5 
has been used. 



1.25(5 — R). While the present paper has dealt with the 
dynamical, or total, mass-to-light ratio (Md yn /Ldi s k) within 
a fixed number of disk scale- lengths (i.e. 4), one can com- 
bine the present observational results with the theoretical 
predictions of Bell & de Jong (2001) to determine how the 

stellar-to-dynamicalpffl| mass ratio must vary with (B — R) 
colour. 

The only strong correlation detected between the dy- 
namical mass-to-light ratio and (B — R) colour is in the 
Tf'-band data of the Cluster Spiral galaxies (Figure 11), 
where log(M dyn /L K , disk ) oc -1.25±0.28(B-i£) (Figure U). 
Combining this result with the stellar mass-to-light ratio 
prediction of Bell & de Jong (2001) (log [M^aiek/Lx, disk] oc 
0.45[B — R]), and neglecting differences between K and K' , 
would imply that the logarithm of the stellar-to-dynamical 
mass ratio must increase as disks become redder, varying as 
1.70±0.28(B — R). If the theoretical models are correct, this 
implies that galaxy modellers should not assume a constant 
stellar-to-total mass ratio ( at least within 4h ) for any given 
passband. This conclusion seems to be at least qualitatively 
correct; LSB galaxies are largely unevolved blue galaxies 
with higher relative gas-to-total mass fractions than HSB 
galaxies (Bothun, Impey & McGaugh 1997). 



II II By stellar-to-dynamical mass ratio, it is meant the ratio of the 
stellar disk mass (excluding any possible bulge mass) to the total 
dynamical mass within 4 scale-lengths. 



7 CONCLUSIONS 

The best-fitting two-dimensional plane to the 3-dimensional 
space of galaxy disk parameters (rotational velocity « ro t, 
disk scale-length h, and central disk surface brightness fio — 
—2.5 log Jo) has been derived for a sample of nearly 100 field 

OC /O.SO±0.05 ft 0.77±0.07 ! and ^ 
0.43±0.04t, 0.69±0.07 



galaxies. In the B-band v Iot ^ j 



the B-band « rot oc ju «±u.u4 /l u.b M ±u.u / Thig implies that in 
the B-band the dynamical M/L ratio is, contrary to pop- 
ular claims, nearly independent of the disk surface bright- 
ness level; the disk scale- length is what is important, with 
log(M/L B ) oc 0.00(±0.04)^ +0.54(±0.14) log ft. This gen- 
eral result holds when using HSB field galaxies, when using 
HSB and LSB field galaxies together, and when derived from 
the TF relation instead of the best-fitting two-dimensional 
plane. We should therefore not talk simply about LSB and 
HSB galaxies, but also give mention to their scale- lengths. 

The M/Lb ratio is shown to be largely independent 
of a disks (B — R) colour, whereas M/L' K varies as - 
1.25±0.28(B — R). Combining this with recent theoreti- 
cal modelling implies that the logarithm of the stellar-to- 
dynamical mass ratio is not a constant value, but increases 
as disks become redder, varying as 1.70±0.28(B — R). Conse- 
quently, galaxy modellers should not treat this as a constant 
quantity. 

In the B-band log(M/L R ) oc 0.06(±0.03)aio + 
0.38(±0.14) log h. Extensions to redder passbands, albeit us- 
ing a sample of Cluster spirals, reveals that the M/L ra- 
tio becomes increasingly less dependent on the disk scale- 
length at longer wavelengths, reaching the situation where 
M/L oc l/V7o at K' . 

Marginal evidence is presented to suggest that the ratio 
of the mass-profile scale-size (o) to the luminous disk scale- 
length (h) may vary with disk luminosity in the B-band, 
such that a/h decreases with increasing disk luminosity. A 
method to test this tentative possibility is described. 
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